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The Nebulosity and Distance of the Cepheid RS Puppis 
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ABSTRACT 

Adopting a distance for RS Pup derived from a period-luminosity relation based on 
Ccphcid parallaxes, it is shown that the phase-lag observations of the surrounding 
nebulosity by Kervella et al. are well fitted by a model of an equatorial disc at an 
angle of 8°.l ± 0°.6 to the plane of the sky. The astrophysical implications of this are 
briefly mentioned. 

Key words: Stars: individual: RS Pup, Stars: circumstellar matter. Stars: distances, 
Stars: variables: Cepheids, ISM: reflection nebulae, ISM: dust, extinct ion 
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1 INTRODUCTION 

Westerlund (1961) discovered that the 41day Cepheid RS 
Pup was surrounded by a remarkable nebulosity. This is in 
the shape of rudimentary rings but with much distorted 
structure and condensations. Havlen (1972) showed that 
portions of the nebulosity varied in the period of the Cepheid 
but with with various phase lags. A very beautiful set of 
measurements of phase lags at various points in the neb- 
ula has recently been obtained by Kervella et al. (2008) {= 
Kervella et al.). In general, the expected phase lag at a point 
i may be written: 



{Ni + A<pi) = 5.775510"^i36'i(l -I- sin ai)/P cos m 



(1) 



Here Acpi is the fractional phase lag, Ni the whole num- 
ber of pulsation periods elapsed, D is the distance to RS Pup 
in parsecs, 6i is the angular distance of i from the star in 
arcsec, P is the pulsation period in days and Qi is the angle 
between the line joining the star to i and the plane of the 
sky (positive if i is further away than the star, negative if it 
is nearer). The measured quantities are A(j>i and 9i. P is as- 
sumed known and here it is taken as 41.4389 days ( Kervella 
et al.). In an attempt to determine D, Kervella et al. assume 
Qi = constant = 0. That is they assume that all the features 
measured by them lie in the plane of the sky and the values 
of Ni are then chosen to obtain the best fit to this model. 
The justification for this assumption is that if the nebulosity 
consisted of a series of thin, uniform, spherical shells centred 
on the star, then the deviation of all measured points from 
the plane of the sky would be small. However an examina- 
tion of the structure of the nebulosity (for instance from the 
figures in Kervella et al.) shows that it is far from corre- 
sponding to this idealized model. There is much distortion 
and density variation in the rudimentary rings. Kervella et 
al. place special emphasis on the ten condensations or blobs 
shown in their fig 7. The existence of such blobs is not consis- 



tent with the idealized model and leaves open the question 
of whether they or other features are actually in, or near, 
the plane of the sky. In view of these uncertainties it can- 
not be claimed that a definitive distance to RS Pup can be 
found based on the "in-the-plane" assumption. In the next 
section this assumption is dropped and it is shown that a 
simple and astrophysically interesting model for the nebu- 
losity is found if a distance for RS Pup is adopted from a 
period-luminosity relation. 



2 AN EQUATORIAL DISC MODEL 

van Leeuwen et al. (2007) established a reddening-free 
period-luminosity relation in V, I based on HST (Benedict 
et al. 2007) and revised Hipparcos parallaxes. This together 
with the data in table Al of that paper leads to a predicted 
distance of 1728pc for RS PupQ Adopting this distance it 
is possible to use eq. 1 to study the three dimensional struc- 
ture of the nebulosity. In principle the values of Ni can be 
arbitrarily assigned. However they should obviously be cho- 
sen to account for apparent continuities in the structure and 
to conform to some simple, physically reasonable model. It 
was quickly found by trial and error that there is a consis- 
tent set of values of Ni in which the points measured by 
Kervella et al. are further away than the star on the south 
side and nearer on the north, i.e. an inclined disc model is 
indicated. This is indeed the simplest model, if the uniform 
spherical shell model is rejected. In such a model the values 
of Ni have to chosen such that {Ni + A(j>i)/9i values are as 
near constant as possible in a given direction from the star 
and vary smoothly with direction. The details are given in 



^ The distance, 1830_g4 pc derived from the pulsational parallax 
by Fouque et al. (2007) is not significantly different from this. 
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Table 1. This contains data on the 31 points observed by 
Kervella et al. and I am greatly indebted to them for sup- 
plying detail of their observations which were not given in 
the original paper. The Table lists: 

1. Position number, i. 

2. Angular distance from the star, di, in arcsec. 

3. Azimuth of the point relative to the star, f3i, measured 
from north through east, in degrees. 

4. A(j>i and its standard error. 

5. The value of Nt (= iV/^) adopted by Kervella et al. to fit 
their model assumptions. 

6. The distance, dy , behind (positive) or in front (negative) 
of the plane of the sky through the star. This is found by us- 
ing eq.l together with adopted values of Ni and D to derive 
at in each case. Then, 
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dy =4.848.10 ^OiDtana^ 



(2) 



where dy is in units of 10^ pc. For dy the value of D esti- 
mated by Kervella et al. (1992 pc) was combined with their 
N^ values. 

7. The value of Ni adopted in the present paper 

8. The distance, dy, behind or in front of the plane of the sky 
through the star assumed to be at its PL distance (1728pc) 
and adopting the revised values of Ni . The units are also in 
10-*pc. 

9. The perpendicular distance dx of the point from the inter- 
section of the disc with the plane of the sky and projected 
onto the plane of the sky. In the same units as dy. This is 
given by: 



dx = 83.776, sin(A - 7) 



(3) 



where 7 is the angle (azimuth) at which the plane of the disc 
cuts the plane of the sky. In this test of the model this is take 
as close to the line from the star to point 9 (i.e. 7 = 80°). 
Fig. 1 shows a plot of dx against dy. This indicates a clear, 
apparently linear relation, between the two quantities. That 
is, the points lie in a tilted plane, presumably an equatorial 
disc. The line shown is a least squares fit through the origin 
and is given by: 



dy = 0.143(±0.010)da: 



(4) 



The tilt of the disc to the plane of the sky is tan^^ 0.143 = 
8°.l ± 0°.6. The rms scatter about the line in dy is only 
73.10~*pc, much smaller than the diameter of the disc, 
which is ~ 6800.10^'*pc out to the limits of the phase- lag 
survey. This rms scatter may be compared to the rms scatter 
of d^ which is lll.lO^^pc. No attempt has been made to 
optimize the disc model by, for instance, varying 7 to find a 
better fit. This might reduce the rms scatter slightly. How- 
ever this is already small compared to the diameter of the 
part of the disc surveyed indicating a relatively thin disc. A 
realistic disc will have some significant depth perpendicular 
to its axis and, indeed, Kervella et al. note that their obser- 
vations at some positions suggest smoothing attributed to a 
non-zero depth in the line of sight. An inclined disc model 
broadly similar the one just discussed could probably be de- 
rived for other distances, if there were good evidence for 
these. However it should be noted that to take the distance 
as a free parameter in an attempt to reduce the rms scat- 
ter in the model is to assume that the disc must conform 
as closely as possible to an idealized model which is per- 



Figure 1. A plot of the distances, dy,from the plane of the sky 
through the star against dx the perpendicular distance, in the 
plane of the sky, to a line of azimuth 7 = 80°. The units are 
lO^^pc. Note the expanded dy scale. 



fectly flat and of negligible thickness. There is no a priori 
justification for such an assumption. 

An equatorial disc model for RS Pup is particularly 
interesting from the astrophysical point of view. Whilst it 
has seemed possible, for instance, that Cepheids might have 
ejected shells in a previous evolutionary phase, it has been 
puzzling that only for RS Pup is such a prominent structure 
found. The interpretation of the nebulosity as a disc at a 
small angle to the plane of the sky opens up the possibility 
of a deeper understanding of this phenomenon and its rarity. 
Obvious possibilities are loss of mass in the equatorial plane 
by unusually rapid rotation and/or binary interaction at an 
earlier evolutionary stage. 



3 CONCLUSION 

The structure seen in the RS Pup nebulosity makes ques- 
tionable the assumption that phase-lag observations all re- 
fer to points close to the plane of the sky and this makes 
distance estimates based on this assumption questionable. 
An inclined disc model at a distance predicted by a period- 
luminosity relation gives a good fit to the data and opens 
new possibilities for understanding the system, including a 
possible interacting binary model. 
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Table 1. The phase lag observations of Kervella et al. with derived linear positions 
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